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The higher order moments of the net-baryon distributions in relativistic heavy ion colhsions are 
useful probes for the QCD critical point and fluctuations. Within a simple model we study the 
colliding energy and centrality dependence of the net-proton distributions in the central rapidity 
region. The model is based on considering the baryon stopping and pair production effects in the 
processes. Based on some physical reasoning, the dependence is parameterized. Predictions for 
the net-proton distributions for Au+Au and Pb-|-Pb collisions at different centralities at ^sjvjv=39 
and 2760 GeV, respectively, are presented from the parameterizations for the model parameters. 
A possible test of our model is proposed from investigating the net-proton distributions in the 
non-central rapidity region for different colliding centralities and energies. 
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I. INTRODUCTION 

The investigation of QCD phase diagram is of crucial 
importance for our understanding of the properties of 
matter with strong interactions. Lattice QCD calcula- 
tions have predicted, at vanishing baryon chemical poten- 
tial, the occurrence of a cross-over from hadronic phase 
to the deconfined quark-gluon plasma phase above a crit- 
ical temperature of about 170-190 MeV [l|,[l. A distinct 
singular feature of the phase diagram is the QCD criti- 
cal point [3] which is located at the end of the transition 
boundary. A characteristic feature of the QCD critical 
point for systems in the thermodynamical limit is the di- 
vergence of the correlation length ^ and extremely large 
critical fluctuations. In ultra-relativistic heavy ion col- 
lisions, however, because of finite size and rapid expan- 
sion of the produced system, those divergence may be 
washed out. As estimated in [s*], the critical correlation 
length in heavy ion collisions is not divergent but only 
about 2-3 fm. So the signals for the critical point of 
the system produced in heavy ion collisions cannot be 
observed as clearly as in the condensed matter physics. 
However, remnants of those critical large fluctuations 
may become accessible in heavy ion collisions through an 
event-by-event analysis of fluctuations in various chan- 
nels of conservative hadron quantum numbers, for exam- 
ple, baryon number, electric charge, and strangeness 
Particularly there would be a non-monotonic behavior of 
non-Gaussian multiplicity fluctuations in an energy scan, 
which would be a clear signature for the existence of a 
critical point. In fact, at vanishing chemical potential it 
has been shown that moments of conservative charge dis- 
tributions are sensitive indicators for the occurrence of a 
transition from hadronic to partonic matter Q. 

Recently, great interest both experimentally Q and 
theoretically J- 9] has been aroused on the higher order 
moments of net-baryon distributions in heavy ion colli- 
sions at the BNL Relativistic Heavy Ion Collider (RHIC) 



energies. The theoretical interest on these higher order 
moments comes from the discovery of the relation be- 
tween the moments and the thermal fluctuations near 
the critical points for the produced quark matter. If some 
memory of the large correlation length in the quark mat- 
ter sate persists in the thermal medium in hadronization 
process, this must be reflected in higher order moments 
of the distributions. Theoretical prediction [ll| showed 
that the third moment, called skewness, is proportional 
to and that the fourth moment, or kurtosis, propor- 
tional to while the second moment proportional to 
More importantly, the moments are closely related to the 
susceptibilities of the thermal medium. Thus the higher 
order moments have stronger dependence on the correla- 
tion length ^ and are therefore more sensitive to the criti- 
cal fluctuations. Recently STAR Collaboration has pub- 
lished experimental data on the higher order moments 
[l^ for different colliding systems at different colliding 
energies for different colliding centralities. It has been 
argued that the net-proton distribution can be a mean- 
ingful observable for the purpose of detecting the critical 
fluctuations of net baryons in heavy ion collisions [ic|. 
This statement makes experimental investigation of net- 
baryon fluctuations much easier, because neutrons and 
strange baryons cannot been detected easily/effectively 
in experiments. Based on theoretical and experimen- 
tal investigations, it has been argued that information 
of QCD phase diagram and the critical point can be ob- 
tained from the energy dependence of those moments fj\ . 

The moments of net-proton distributions have been 
studied recently by quite a few groups with differ- 
ent event generators such as A Multi-Phase Transport 
(AMPT) and Ultrarelativistic Quantum Molecular Dy- 
namics (UrQMD) [3, Heavy Ion Jet INteraction Gen- 
erator (HIJING) and hadron resonance gas model 
[m, 15] etc. Some other authors tried to search the sta- 
tistical and dynamical components in the net-proton dis- 
tributions in [l6| . where the statistical distributions for 
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both proton and anti-proton are assumed Poissonian and 
the departure from Poisson distributions is regarded as 
the dynamical influence. One should pay attention that 
the use of independent Poisson distributions for proton 
and anti-proton implies that protons and anti-protons are 
produced completely uncorrelated. Therefore the baryon 
number conservation may be violated in any event. On 
the multiplicity distribution of hadrons, a canonical en- 
semble is employed in jl7j to derive the number distri- 
bution for TT systems. This is reasonable because there 
are a lot of tt particles in the final state of heavy ion 
collisions. But a simple transportation of the method to 
the case for baryon production may be problematic, be- 
cause the relevant baryon particle number may be not 
large enough for an equilibrium statistical description. 
In Ref. [ISi] the net-proton distributions in Au-|-Au col- 
lisions at ^/SNN — 200GeV are studied from very sim- 
ple but well established physics considerations: baryon 
stopping and baryon pair production. For a given mean 
net-proton number from initial nuclear stopping, the ini- 
tial proton number is assumed to satisfy a Poisson dis- 
tribution. From the produced baryon pairs, the joint 
distribution for the newly produced proton and anti- 
proton can be derived. Then one can obtain the net- 
proton distribution in the final state of heavy ion colli- 
sions. Good agreement with experimental data has been 
obtained for Au-|-Au collisions at three colliding central- 
ities at ^/snn — 200GeV. 



This paper is an extension of the work in Ref. 18 1 
by studying the moments of net-proton distributions in 
a given central rapidity window in Au-|-Au collisions at 
lower RHIC energies at different centralities. This paper 
is organized as follows. In next section, we will address 
our model and the physics points for the centrality de- 
pendence of parameters. Using an analytical expressions 
for the net-proton distribution derived in Ref. |18j we will 
show that the moments of the distribution up to fourth 
order can all be well described by our model with suitably 
chosen parameters for four colliding energies in Au-|-Au 
collisions at RHIC. In section III, our model results for 
the moments are compared with the experimental data 
from STAR Collaboration. The net-proton distributions 
for Au-|-Au collisions at ^snn=39 GeV are presented 
for three centralities. Then in section IV, we discuss 
the energy dependence of the parameters and give pa- 
rameterizations for such dependence. Then we extrapo- 
late the dependence to the CERN Large Hadron Collider 
(LHC) energy and predict the net-proton distributions 
for Pb+Pb collisions at different colliding centralities for 
y^SNN =2.76 TeV. The last section will be for a brief 
summary. 



II. MODEL CONSIDERATION FOR THE 
CENTRALITY DEPENDENCE 

Nuclear stopping plays an important role in heavy ion 
collisions and the study of such effect is a fundamental 



issue, since this effect is related to the amount of energy 
and baryon number that get transferred from the beam 
nucleons into the reaction zone. We denote B the mean 
net-proton number in the final state distribution. As can 
be seen from our model consideration, B comes only from 
the initially stopped protons. In nuclear-nuclear colli- 
sions the mean net-baryon number B in central rapidity 
region would be zero if there were no nuclear stopping 
in the processes. Another physics point we consider in 
Ref. |18j and here is baryon pair production in the in- 
teractions. We denote ^ the mean number of produced 
baryon pairs within a given kinematic region in the col- 
lisions at given colliding centrality. By assuming that 
baryon pairs are produced independently, the pair num- 
ber distribution is of Poissonian. With isospin conserva- 
tion, one can derive a simple analytical formula for the 
net-proton number Ap distribution P(Ap) as a function 
of B and fj, as [l3| 



P(Ap) 



^ dx 



Now we study the centrality dependence of B and /i. 
Because the nuclear stopping effect results from the inter- 
actions between a passing nucleon with other nucleons on 
its way in a nuclear-nuclear collision, the baryon number 
stopped in a given kinematic region is closely related to 
the number of participant nucleons Np. In more central 
collisions the stopped net baryon number will be larger. 
In nuclear- nuclear collisions, if every nucleon from a nu- 
cleus suffers exactly the same interactions, the stopped 
proton number would be -B oc iVp. Of course, the real 
case is not so simple. Because multiple scattering ef- 
fect is more important for more central collisions, a little 
larger B/Np can be expected for more central collisions. 
On the other hand, the stopped net-proton can be de- 
tected only in the final state, or in other words, after 
evolving with the system for some time. During the evo- 
lution of the system, the net-proton may diffuse into a 
kinematic region out of our interested window. For cen- 
tral collisions, the evolution time is longer and such dif- 
fusion effect is more obvious. This effect would reduce 
B/Np for central collisions. The real proportional factor 
B/Np is a result from the competition of the two effects: 
initial multiple scattering and later baryon number dif- 
fusing. Therefore, the factor B/Np should have a weak 
centrality dependence. Thus one may parameterize the 
centrality dependence of B as 



B = aiNp{l - aaiVp) , 



(2) 



with ai and a2 depending on the colliding energy of the 
system. While ai is always positive, the magnitude of 
a2 should be small and a2 can be negative or positive, 
depending on whether multiple scattering is more im- 
portant than baryon number diffusion or the opposite. 
For the baryon pair production, similar physics consid- 
erations apply also. The probability for a baryon pair 
production near a point in the system is determined by 
the energy density at that point. If the energy density 
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V«NN (GeV) 


ai 


lO^'aa 


&i 


lO-^ba 


19.6 


0.022 


-0.126 


0.0132 


0.121 


39 


0.0012 


-1.97 


0.0300 


4.90 


62.4 


0.0096 


-0.159 


0.0383 


5.10 


200 


0.0052 


1.88 


0.0585 


5.00 



TABLE I. Fitted parameters for four colliding energies. 



of the hot strong interacting matter is uniform and the 
same for all colliding centralities, one may expect pro- 
portional to the volume of the system, thus ji cx iVp. But 
the initial energy density is higher for more central col- 
lisions due to stronger nuclear stopping effect, a larger 
fJ./Np is expected for more central collisions, if the nu- 
clear stopping is the only physics in the process. On the 
other hand, as in the consideration for i?, the energy 
diffusion from central to non-central rapidity region will 
reduce the value of fi/Np. The competition of these two 
effect results in a behavior of /z as a function of A'p simi- 
lar to that of B. Therefore we parameterize the colliding 
centrality dependence of fj, as 



H = hNp{l - b2Np) 



(3) 



with bi and 62 also depending on the colliding energy. 
As for 02, the value of 62 should be very small but can 
be positive or negative due to the competition of initial 
energy stopping and diffusion in the evolution of the sys- 
tem. 

With the above four parameters, one can calculate the 
net-proton distribution and all the associate moments for 
any colliding centrality for given center of mass energy of 
the colliding system. 



III. COMPARISON WITH THE 
EXPERIMENTAL DATA 

For Au-f An collisions at RHIC energies, we investi- 
gate the moments up to fourth order for the distribution 
of net-proton in the central rapidity window \y\ < 0.5 
as functions of iVp by using our model described in the 
last section. The fitted results from our model for the 
moments are shown in Figs. 1-4 for the mean, vari- 
ance, skewness and kurtosis for four colliding energies 
ysii7=19.6, 39, 62.4 and 200 GeV. The fitted parameters 
are tabulated in TABLE L It can be seen that our simple 
model can describe quite well the centrality dependence 
of moments for the four energies with the parameters 
chosen. From the excellent agreement with the experi- 
mental data, one can conclude that our model contains 
the necessary physics for the net-proton distributions. 

With those parameters in TABLE HI one can calcu- 
late the net-proton distributions quite easily for differ- 
ent centralities for the four energies as in TABLE HI Our 
new parametrization for B and can give distributions 
for Au-t-Au collisions at ■\/snn=200 GeV. The newly ob- 
tained distributions have no visible difference from those 
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FIG. 1. Mean value of the net-proton distributions as a 
function of A^'p for Au-|-Au collisions at four ^snn, (a) 200; 
(b) 62.4; (c) 39 and (d) 19.6 GeV. The points are from 
RHIC/STAR data [H. 
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FIG. 2. variance of the net-proton distributions as a function 
of Np for Au-|-Au collisions at four y^SNN, (a) 200; (b) 62.4; 
(c) 39 and (d) 19.6 GeV. The points are from RHIC/STAR 
data [11 . 



in Ref. [l8||, in good agreement with the STAR data, so 
will not be presented here. As an example to show the 
distributions, we present here only the distributions at 
Y^SNN = 39 GeV for Au-I-Au collisions at different collid- 
ing centralities, for later comparison with experimental 
results. 



IV. COLLIDING ENERGY DEPENDENCE OF 
NET-PROTON DISTRIBUTION 

After discussing the net-proton distributions for 
Au-|-Au collisions at four different colliding energies, one 
can discuss the colliding energy dependence of parame- 
ters in our model. As we discussed in Sec. II, the values 
and their dependence on the colliding energy can tell us 
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FIG. 3. Skewness of the net-proton distributions as a function 
of A^p for Au+Au collisions at four ^snn, (a) 200; (b) 62.4; 
(c) 39 and (d) 19.6 GeV. The points are from RHIC/STAR 
data [11. 
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FIG. 4. Kurtosis of the net-proton distributions as a function 
of A'p for Au+Au collisions at four ^/sm^, (a) 200; (b) 62.4; 
(c) 39 and (d) 19.6 GeV. The points are from RHIC/STAR 
data El. 



some physics in the colliding process, especially the com- 
petition of effects from the initial multiple scattering and 
later baryon number transportation. With the increase 
of colliding energy, the stopped baryon number in the 
central rapidity region will decrease. Thus the value of 
ai will decrease with y/swN- Though the energy fraction 
stopped in the central rapidity region becomes smaller 
at higher colliding energy, the energy density in that re- 
gion still increases with ^/sjm- Then bi will increase with 
Y^snn- For y/swN high enough, one may expect ai and 
bi saturates at some limiting values. The behaviors of 
a2 and &2 can be quite different from those of ai and 
bi. Since 02 and 62 depend on the competition of effects 
from initial multiple nucleon-nucleon scattering and later 
baryon/energy diffusion, the behaviors of 02 and &2 with 



FIG. 5. Net-proton distributions for Au-|-Au collisions at 
y^SNN = 39 GeV at different colhding centralities. 



the increase of ^snn can be complicated. Both the ef- 
fects from the initial multiple nucleon-nucleon scattering 
and later baryon/energy diffusion become stronger with 
the increase of y^snn, but their rates of increase may be 
different. For a given increase of ^snn, if the initial mul- 
tiple scattering becomes more important, |a2| (I62I) will 
increase with -^snn in the region with negative 02 (&2)- 
Otherwise, \a2\ (I&2I) will decrease in the same region. 
For y/swN — > 00, the interaction duration in the produced 
matter will very long and the diffusion effect will be much 
stronger than that from initial multiple scattering. Then 
one can expect 02 and &2 approaching some positive sat- 
urating values, implying that B/Np and ^i/Np is smaller 
for central collisions when ^snn is high enough. 

To see the colliding energy dependence of the four pa- 
rameters, we plot the fitted parameters shown in TABLE 
I as functions of y^snn in GeV. The plots are shown in 
Figs. [n]and[71 The shown dependence of the parameters 
on the colliding energy can be described by the following 
expressions 



10^02 

bi 
1046, 



14.55(1 + 9.23 X 10-3^/ii;^)/(l + 40.2yiii) , 
1.87-1.21 X 10-''(7iii7)'^exp(-0.1077iii') 
-0.02(1 - 0.10575ii7)/(l + 0.029575ii7) , 
5.0- 0.626(7iii)^ "exp(-0.9797si^) . 

(4) 

The functional form for ai and bi {i — 1,2) are chosen to 
satisfy the demands from the physics considerations in 
the last paragraph. With the above expressions for the 
energy dependence of the parameters, it is straightfor- 
ward to calculate the values of those parameter at the 
LHC energy ■y/SNN=2760 GeV. By assuming that our 
model can be applied to that energy, one can calculate the 
net-proton distributions at that energy for different col- 
liding centralities. The obtained distributions are shown 
in Fig. [S] The values of Np used in the calculation are 
from p|. 

In all considerations up to now, we only discussed the 
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FIG. 6. Colliding energy dependence of parameters ai and 02 
for the initially stopped proton number in the given central 
rapidity window. Points marked by star are from our model 
fitting, and the points marked by triangle are calculated from 
Eq. (4) at ^/s^=2760 GeV. 
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FIG. 7. Colliding energy dependence of parameters bi and 62 
for the mean number of produced baryon pairs in the given 
central rapidity window. Points marked by star are from our 
model fitting, and the points marked by triangle are calcu- 
lated from Eq. (4) at y/s^=2760 GeV. 



net-proton distributions in the central rapidity region 
\y\ < 0.5. Because of the diffusion of baryon number 
and energy from central to non-central rapidity region, 
the parameters 02 and 62 show complicated behaviors 
as functions of colliding energy. If we consider the net- 
proton distributions in the non-central region, \y\ > 0.5 
for say, the same physics arguments apply and one can 
expect that the distributions can also be described by 
Eq. (1) with centrality dependence of parameters B and 
H being given by Eqs. (2) and (3). One can also expect 
that the colliding energy dependence of ai and bi is sim- 
ilar to that for the case in central rapidity region. For 
a2 and 62, the baryon number and energy diffusion from 



central to non-central r(\2;ion will not compete to but co- 
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FIG. 8. Net-proton distributions for Pb-f-Pb collisions at 
^/snn = 2760 GeV at different colliding centralities. 

operate with the multiple scattering effect. Then 02 and 
&2 will be negative for all colliding energies and all A-l-A 
collisions. This statement can be tested experimentally. 

V. CONCLUSION 

We studied the net-proton distributions for Au-t-Au 
collisions at four colliding energies for ^snn from 19.6 to 
200 GeV at different centralities. Based on some physical 
arguments, the parameters in our model are parameter- 
ized as functions of centrality and energy. The higher 
order moments for the distributions are in good agree- 
ment with the experimental data. Prediction for the net- 
proton distributions at LHC energy y^snn=2.76 TeV is 
presented for different centralities. The net-proton dis- 
tribution in a non-central rapidity region and its depen- 
dence on centrality are discussed. 

It should be mentioned that nothing else is assumed 
in this model except an initial stopped net-proton and 
a finite probability for producing baryon pairs from the 
produced matter. Therefore, our model has nothing 
to do with thermal equilibrium and/or critical fluctua- 
tions. Because our model consideration is based on nor- 
mal physics effects, our results can be used as a baseline 
for detecting novel physics in the processes. 
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